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an amazing history, and it might influence our perception of a symmetry 
beyond this example. 
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1 Introduction 



After the introduction of non-Abelian gauge theories ^ it took a long time to 
understand how to set up perturbation theory |2] and control the ultraviolet 
fluctuations [3| in such models. Asymptotic freedom 0] raised the possibility to 
connect non-Abelian gauge theories with the phenomenology of strong interac- 
tions. It was not clear, however, what to do with the wild infrared divergences 
present in perturbation theory. A few courageous theorists with good intuition 
interpreted this problem as a virtue: objects with color are confined by strong 
forces and the infrared problem in perturbation theory is an indication of that. 
Quantum chromodynamics (QCD), the theory of quarks in interaction with glu- 
ons in an SU{3) gauge theory became an attractive candidate to describe the 
physics of hadrons ^ . 

In 1974 Wilson W presented a,U{l) gauge theory with fermions on a hyper- 
cubic four dimensional Euclidean lattice which has confinement in the g —^ oo 
limit, where g is the bare gauge coupling. Although in this limit the model lives 
on a coarse lattice far from the continuum this work introduced many of the 
ideas, tools and notions which lead to a new field in particle physics. The formu- 
lation is fully gauge invariant and independent of perturbation theory. It creates 
also a strong link to critical phenomena in statistical physics. The construction 
was generalized to SU{N) in the same paper. The confinement problem can 
be formulated in this language the following way: do these models stay in the 
confining phase as the bare gauge coupling is tuned towards g = where the 
continuum theory is defined? It took a few years until the new tools (mean field 
and strong coupling expansion techniques [7|, Hamilton formulation |S], Monte 
Carlo simulations jHI) have been adapted to quantum field theories, in particular 
to QCD. 

In 1980 Creutz demonstrated ^01 in a Monte Carlo simulation that close 
to the continuum limit the string tension a between a static quark-antiquark 
pair in an SU{2) Yang-Mills theory behaves as expected from continuum renor- 
malization group considerations. This was a spectacular result indicating that 
one can get close to the continuum with modest computing power in this non- 
Abelian gauge theory and confinement survives this limit. Of course, today's 
state-of-the-art calculations control the systematic and statistical errors of the 
static quark-antiquark potential both in SU(2) and SU{3) gauge theories much 
better. The potential is followed up to 4fm distances, the fine details of the fluc- 
tuating color flux tube are seen and quantitatively determined 11 . Although it 
is not proven rigorously, there remains little doubt that in these gauge theories 
the static quarks are confined in the continuum limit. 

QCD, the theory of strong interactions is, however, a Yang-Mills gauge the- 
ory coupled to dynamical quarks. It has an extremely rich phenomenology most 
of which is beyond the reach of perturbation theory. QCD is an exciting theory 



in itself. In addition it enters most of the processes in the standard model in 
a relevant way either in hadronic matrix elements, or as a background. For 
a recent general overview on the theoretical and numerical activities in lattice 
QCD I refer to ref. ^21 and the references therein. 

In this contribution I will mainly consider the history of chiral symmetry 
on the lattice which was long and troubled and had an amazingly nice upshot. 
This outcome might influence our perception on the realization of a symmetry 
in a quantum field theory in general. The following overview is non-technical. 
For further details and applications I refer to the summaries in |13[ I14L [T5| . 



2 QCD on the lattice 



A hyper-cubic lattice is a natural regularization when path integrals are used 
to formulate a quantum field theory in d = 4 Euclidean space. The quark 
field ip lives on the lattice points indexed by an integer vector x G Z4 while 
the gauge fields are associated with the links (x, /i), /.t = 1, . . . , 4 of the lattice. 
Derivatives are replaced by finite difference operators using, for example the 
nearest neighbor forward and backward difference operators^ 

^K.a ~ "x+il.y — ^x,y, C)x,y ~ "x,y ~ "x-p.,y ■ (1) 

A simple guess for the action of a massless free fermion on the lattice might 
read as 

x,y 
Switching on the gauge interaction a parallel transporter is included in the finite 
difference operators to assure gauge invariance 

Vx,y = U{X, H) 5x + [i,y - 5x^y , V'^*y = S^^y - U{X - fl, fiY S^^f^^y , (3) 

where U{x, fi) is an element of the gauge group SU{3). The gauge field U{x, fj) 
can be expressed in terms of the vector potential A^ [x) used in the continuum 
formulation as 

U{x,ii) = Pexp[i I drA^'ix + (f - t)(i)] , (4) 

Jo 

where P denotes r-ordering and A^ is an element of the Lie algebra. On the 
lattice, however, the group element U{x,fj,) is the fundamental field variable 
which transforms under a gauge transformation as 

U{x,fi)^V{x)U{x,i2)V\x + fi). (5) 



■^Here and in most of the following expressions we suppress the lattice unit a. All the fields 
and parameters are made dimensionless by appropriate a powers: in the mass m a factor of 
a, in the fermion field ^ a factor of a^'^, etc, is absorbed. 



Eq. lO implies that the trace of the product of gauge matrices along a closed 
path on the lattice is a gauge invariant quantity. Locality requires that the 
extension of the gauge action is 0{a). Taking the smallest loop, the plaquette, 
we get a discretized, local and gauge invariant action for QCD 

'^ = -4E^«Trf/p + ^^,i?,,,Vy, (6) 

P x,y 

where the sum runs over the plaquettes, Up is the product of four directed link 
matrices around the plaquette p and the Dirac operator has the form 

^= ^T^IV^ + V^*) . (7) 

Indeed, taking the formal continuum limit a — > of eq. © and using eqs. (|3I4|I 
we obtain the standard continuum form of the QCD action for a massless quark. 
Obviously, there are infinitely many local and gauge invariant lattice actions 
which have this property. The continuum predictions must not depend on this 
freedom, however. In an asymptotically free theory the scaling dimension of 
operators agrees with their engineering dimension up to logarithms (at least in 
perturbation theory) and it is easy to count the number of relevant operators. 
Perturbative renormalizability can be demonstrated on the lattice also .16) . 

In order to define the quantum theory we have to fix the integration measure 
of the path integral. Here we take the natural choice of the group invariant mea- 
sure for every link variable U{x,y) and integrate over the Grassmann variables 
^(a;) and 'ipix) according to the standard rules of Grassmann integration. This 
measure is also gauge invariant. 

The action in eq. IJBJ has a problem, however. The Dirac operator in eq. lO 
describes actually 16 massless fermion species - a problem called doubling. This 
is the first small surprise along a tortured way around chiral symmetry. 



3 Doublers, chiral symmetry and a 
no-go theorem 

The free fermion Dirac operator D in Fourier space has the form 

which has 16 zeros in the first Brillouin zone at (0, 0, 0, 0), ... , (tt, tt, tt, tt). It can 
be shown |E1 that these species couple to the axial current with alternating sign 
and their total contribution to the U{\) axial anomaly cancels |18lll9j (see also 



the contributions of Adler and Jackiw in this volume). This observation resolves 
a paradox: the theory is fully regularized, D anticommutes with a standardly 
defined 75, so we have a U(l) chiral symmetry without anomaly. Without the 
doublers this would be a contradiction. 

The question remained whether there exists a clever lattice discretization 
such that the doublers are removed but (at least the non-singlet part of) chiral 
symmetry survives on the lattice. This question was answered by a no-go theo- 
rem l^ni (see also ref. ^7]). The no-go theorem states that for free fermions the 
following four conditions can not hold simultaneously: 

1. D{x) is local; 

2. the Fourier transformed D behaves for p ^ 1 as 17''^^ + 0{p'^) ; 

3. there are no doublers; 

4. 75D + ^75 =0 . 

If we do not want to violate flavor symmetry for massless fermions then 
D{x) is diagonal in the flavor indices and the statement is valid for each flavor 
separately. Locality is the most important property of a quantum field theory, 
so we can not give up this condition. The 2nd point just declares that we have 
a Dirac particle. If we do not want doublers (since there are no such copies in 
nature), there remains only the possibility to violate condition 4, i.e. to give up 
chiral symmetry on the lattice. 

Accepting chiral symmetry breaking terms in the action there are different 
possibilities to follow. Wilson [2] suggested to add a dimension 5 operator to 
the action in eq. JHJ which modifies the Dirac operator as 

£'w = ^[7^(V^ + V^*)-V^V^*] . (9) 

The new term (Wilson term) gives large oc cutoff masses to the doublers, but 
- at least classically - leaves the p — Q pole unchanged. The Wilson term has 
an effect on the renormalization of bare parameters (in particular, for the bare 
fermion mass niq = is not protected), creates mixings between operators in 
different chiral multiplets and it also influences the way and speed the continuum 
is approached. The final physical predictions of the continuum theory are, 
however, independent of this term. 

One might keep a part of the chiral symmetry by giving up fiavor symmetry 
on the lattice [SlEl]. The Kogut-Susskind (or 'staggered') fermions describe 4 
fiavors and keep a (non-singlet) U(l) part of the original C/(4) x C/(4) symmetry. 
This formulation has many attractive properties, in particular for numerical 
simulations 23^. On the other hand, getting QCD with 3 light flavors is non- 
trivial with staggered fermions. I shall not follow this interesting possibility 
further in this contribution since staggered fermions are not directly related to 
the developments I want to discuss. 



4 The Ginsparg-Wilson relation 



The no-go theorem seemed to close all the ways towards a chiral symmetric 
lattice regularization. As the authors formulated ^^: 'The important conse- 
quence of our work is to discourage any attempt to construct chiral invariant 
lattice models for QCD\ No-go theorems, however, are frequently circumvented 
in an unexpected way. This is what happened with chiral symmetry on the 
lattice. 

Soon after the no-go theorem was presented Ginsparg and Wilson 24, sug- 
gested a condition for the Dirac operator which, as the authors formulated, 
implies a 'remnant' chiral symmetry on the lattice. More than 15 years later 
this condition, called the Ginsparg-Wilson relation, became very important in 
the theoretical developments concerning global and local chiral symmetry on 
the lattice and beyond. 

Ginsparg and Wilson treated free fermions and used some basic notions 
of Wilson's renormalization group theory. Consider a renormalization group 
transformation which blocks the fcrmion field ipx living on the fine lattice into 
the Xxb fermion field on the coarse lattice: 

eM-Ssix, X)) = J Dl^D^ exp[ -5(V^, V) 



xb,Vb 



(10) 



Here S and Sb are the original and blocked actions, respectively, uJx,xb describes 
the averaging process in a local neighborhood of the coarse lattice point xs, 
while R is an arbitrary local operator whose inverse is also local and diagonal 
in Dirac space''. Take the fine lattice infinitely fine, so the field -0 lives in the 
continuum and S" is a chiral invariant action in the continuum. Since the block 
transformation explicitly breaks chiral symmetry (it has &X---X structure), so 
does the action Sb on the lattice. The renormalization group transformation, 
however, does not change the long distance (in lattice units) behavior: Sb 
should remember that the starting action was chiral invariant. It is a few steps 
of algebra to dig out this information |23 . Writing 

Sb= ^ Xxb DxB,yB Xys (H) 

XB,yB 

the Dirac operator is constrained by the relation 

Dj^ +-/'^D = Dj^2RD. (12) 

•^The simplest example is Rxb-vb ~ 2^^B,yB 



The physical imphcations of the Ginsparg- Wilson relation in eq. 1121) which 
was obtained for free fermions go far beyond of what the derivation above might 
suggest. Assume there exists a local, gauge covariant solution of this equation 
in lattice QCD. The first trivial observation is that the inverse of D satisfies 

i:i-i-y5+-^5^-i =752i?. (13) 

The inverse of D, the quark propagator, from which the physical correlators are 
constructed is a non-local quantity. The equation above shows, however, that 
its anticommutator with 7^ is local, i.e. has an extension of 0{a) only. The 
violation of chiral symmetry is a contact term. It has been shown in ref. |23 also 
that a Dirac operator satisfying eq. H12() reproduces the U{1) chiral anomaly. 

The following simple consideration strongly suggests that lattice QCD with 
a Dirac operator satisfying eq. (|12|l gives chiral invariant answers for all the 
physical questions. Consider an arbitrary hadron correlator where the opera- 
tors are separated by physical distances. Write down the corresponding Ward 
identity using the standard method in path integral formulation: introduce new 
Grassmann integration variables which are related to the original ones by an 
infinitesimal 7^ chiral rotation. The action is not invariant and, according to 
eq. (|12|l . gives an insertion oc xD^^2RDx to the correlation function. The 
fermion field x (x) will be paired with one of the antifermion (fermion) fields 
of the hadron operators. The two propagators which bind the breaking term 
to the hadron operators are canceled by the two D operators in the breaking 
term. What remains is ^^2Rx,y, where x and y are the positions of the hadron 
operators. Since R is local and the hadron operators are separated at physical 
distances, the contribution from the breaking term is zero. 

Eq. H12|l is a highly non-trivial equation. The chance of finding a local 
solution in the interactive case was rather dim. The paper and its results were 
hardly noticed and got completely forgotten^. 



5 Improvement programs 



On a lattice with finite (in physical units) lattice spacing a the predictions 
systematically deviate from the final a — > continuum values. To increase the 
resolution (decrease the cutoff effects) by a factor of two in a given physical 
volume requires at least 16 times more computing work and memory. Since 
the final continuum limit is universal but the cutoff effects are not, the search 
for lattice actions with reduced cutoflF effects is an important issue since the 
beginnings. Although this question was raised from the numerical side, the 

^According to SLAC Spires the paper was not cited at all over twelve years between 1986 
and 1997 September when the first solution of the GW equation in the presence of gauge fields 
was identified. 



theory behind the improvement programs reveal interesting aspects of quantum 
field theories and it is also part of the developments around chiral symmetry. 



5.1 Symanzik improvement 

To illustrate the basic idea consider a free massless scalar field with a simple 
nearest neighbor action on the lattice 

^0 = E \dt.^{^)d^.cl>{x) . (14) 

X 

The pole of the propagator in Fourier space defines the energy as the function 
of the momentum. For small momenta the energy has the form 

E^=v'-y^Uf+JZvt\+0{a^), (15) 

where for better visual understanding the lattice unit a is written out explicitly. 
The leading cutoff effect is Ola?), ft is easy to see that by adding a dimension 
6 ('irrelevant') term cSi to 5o, where 

Si^Yl 2^f^9f,(l3{x)df,d^(l){x) (16) 

and ci = 1/12 the 0{a'^) cutoff effect cancels. The action S = Sq + ciSi is an 
0{a^) improved lattice action which approaches continuum limit significantly 
faster than the original action Sq. 

In the interacting case the situation is more involved. As it is well known 
the coefficients of the higher dimensional operators enter the relation between 
the renormalized and the bare parameters. This relation contains diverging 
powers and logs of the cutoff. Expressed in terms of the renormalized parame- 
ters, however, the continuum predictions show no sign of the presence of higher 
dimensional operators (universality). They have a role, however, in the improve- 
ment. 

In spite of the technical complications, the final statement on improvement 
is bold and simple: by adding to the original interactive action a linear combi- 
nation of the dimension 6 operators with properly chosen coefficients, the 0{a^) 
cutoff effects can be eliminated in all physical quantities to all orders of pertur- 
bation theory. This has been shown by Symanzik [21] for the (/>'* scalar model 
and in the asymptotically free d — 2 non-linear a-model. This is a highly non- 
trivial result which has been obtained with the help of Callan-Symanzik type of 
renormalization group equations and local effective lagrangians. 



The improvement technique can be extended to gauge theories and QCD |26| . 
During the last two decades large effort was invested to calculate the improve- 
ment coefficients not only for the action but for currents and other densities 
also. These coefficients can be calculated in perturbation theory, or by simula- 
tions non-perturbatively. The Symanzik improvement has an important role in 
controlling the cutoff effects in different applications [27| ■ 

In QCD, if chiral symmetry is broken on the lattice, the leading cutoff ef- 
fects are 0(a) which makes improvement even more important. 0{a) improve- 
ment can be achieved by adding a lattice version of the dimension 5 operator 
ipcr^i,F^i,ip with a tuned coefficient |26j . This term breaks chiral symmetry and 
has no place in a chiral symmetric formulation. A chiral symmetric action is 
automatically 0{a) improved making this possibility even more attractive. 

The C/(l) (non-singlet) chiral symmetry preserved by the 4-flavor staggered 
fermions assures 0(a) improvement. Nevertheless, as painfully realized, this ac- 
tion has large cutoff effects on lattices which are typically simulated today |28| . 
Although the problem could be alleviated by canceling Oi^a^) flavor symme- 
try breaking effects ^^, the message remains. Close to the continuum limit 
0{a) {0{a'^)) cutoff effects dominate in bosonic (fermionic) theories. Whether 
in today's simulations this is the case depends on the theory, the space-time 
dimension, the form of the leading action and on the form of the higher dimen- 
sional operators chosen. Accidentally, just the d = 2 asymptotically free 0(3) 
non-linear cr-model, which was used by Symanzik to illustrate how to extend 
his method to theories with constraints, defies compliance with the Symanzik 
program in numerical simulations |H()j . 



5.2 Classically perfect actions 

There exist lattice actions which have no cutoff effects whatsoever in the classical 
theory. Take the d ^ 2 asymptotically free 0(3) non-linear (T-model mentioned 
above. One can put this theory on a quadratic lattice with a local, classically 
perfect action. The corresponding Euler-Lagrange equations have exact scale 
invariant instanton solutions, they carry a topological charge Q which is an 
integer, on any configuration the value of the action S > 47r|Q| and S = 47r|(5| 
on solutions. If the interaction is switched off the particles have the exact 
continuum dispersion relation E — |p|. The existence of such an action follows 
from Wilson's renormalization group (RG) theory: the classically perfect action 
is the fixed point (FP) of a RG transformation and is determined by classical 
saddle point equations ^^. 

The same is true for QCD in d = 4 [SI ESI- Write the QCD action in the 
form (3Sg + Sj, where Sg and Sj denote the gauge and fermion parts, /3 oc l/g"^ 



and consider the RG transformation in QCD 
exp[-{(3'S'g{V) + S'^{x,X,V))] = 

Di;Di;DUe^p[-p{Sg{U)+Tg{V,U)) + {Sf{i>,i^,U)+Tfix,X;i^,i^,U))]. 

(17) 

Here Tg{V, U) defines the gauge invariant averaging process from the fine U to 
the coarse field V . The averaging function for the fermions Tf{x, X! "0: V': U) can 
be chosen in the form like in eq. H1U|) after introducing parallel transporters when 
adding fermion fields in different lattice points in the averaging: oj — > oj{U) and 
similarly R^^ -^ R^^{U). There is a considerable freedom here, the details are 
not important. In an asymptotically free theory the FP is at /? ^ cxd where the 
path integral above is reduced to classical saddle point equations. For the gauge 
action this reads as 

5f (F) = min[5f (f/) + Tg{V, U)] , (18) 

while for the fixed point Dirac operator in S^^ one obtains 

^^''(^)xi,,« = RiV).,.ys + E ^(C/).B.xZ?^P(t/)-XC/)t ,^ . (19) 

n,n' 

Using the FP equation eq. ((T5)l for the gauge action S'J^ one can derive re- 
sults similar to that discussed above for the 0(3) a-model on classical solutions, 
scale invariant instantons, etc. at finite lattice resolution. The FP Dirac oper- 
ator from eq. (|19|l has exact dispersion relation in the free field limit and gives 
the correct continuum magnetic moment (g-factor) independently of the lattice 
resolution. 

The action S^^ is defined if for any configuration V the value S^^{V) (a 
real number) can be calculated. Take a configuration V on an L* lattice and 
assume that the blocking is from L' = 2^L to L. Considering any valid (gauge 
invariant, local) gauge action on the r.h.s. of eq. (|18|) the minimum gives the 
value of 5'J^(F) with an error which goes to zero as ex A~^ . If we insert an 
approximate solution for the FP gauge action on the r.h.s., the error is reduced 
further. The convergence of the solution of the FP Dirac operator in eq. p9|l is 
oc2-'=. 

Considering the action (38^^ + ^D^^^ for finite /3 we are off the renormal- 
ized trajectory (where the action would be quantum perfect). It is expected 
nevertheless that a classically perfect action will perform well in quantum sim- 
ulations also. This is demonstrated in different d = 2 and 4 models including 
SU{3) Yang-Mills theory and in (quenched) QCD pi|. 

In a numerical simulation, where we need the value the gauge action and 
the D^^v matrix-vector product frequently, there is no way to calculate them 
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using the FP equations above. One can, however, parametrize the solutions in 
terms of gauge loops for Sg and a certain number of fermion offsets and paths 
for D, fit the coefficients to eqs. (|18ll9f) and obtain an approximate solution. 

This improvement program is less systematic than Symanzik improvement. 
Very close to the continuum limit the Symanzik program with non-perturbatively 
determined coefScients is the best procedure. On today's lattices, however, the 
approximate FP actions in many cases perform better than Symanzik improved 
actions |23]. They are, however, also more expensive. 

As I shall discuss, the FP Dirac operator played a relevant role in the quest 
for a chiral symmetric lattice regularizations also. 



6 Domain-wall and overlap fermions 



In 1992 Kaplan |25I suggested a new approach for lattice fermions allowing them 
to move in a five-dimensional space. Setting up a four-dimensional domain- wall 
some light modes become localized to the wall, the other, non-localized modes 
remain heavy on the cutoff level. The idea is related to earlier considerations in 
the continuum |36| . 

Consider a d = 5 free Dirac operator in the continuum with an 2:5 dependent 
mass term 

D = Yd^ + "f^d"^ - M{xr,) , (20) 

where M{x^ -^ ±00) — ±n^o and forms a kink between these values at ^5 = 0. 
The Dirac equation has a massless (ip4 = |p|,p = (pi,P2,P3)) left handed 
solution which is localised to the d = 4 kink at X5 = 0. All the other modes 
live on the cutoff level ttiq. Kaplan realized that this construction can be taken 
over to the lattice and the mechanism works in the presence of gauge fields also. 
Actually, the five-dimensional model is not a genuine gauge-fermion theory: 
the gauge fields know nothing about the fifth direction. It might be better to 
consider this extra dimension as an internal 'flavor' space and the associated 
new fermion fields as regulators which are there to keep chiral symmetry |37) . 
In this approach, which is related to other descriptions [32] also, the effect of 
the massless chiral fermion can be represented as an overlap of two fermionic 
states P7] . 

The domain-wall approach for vector theories like QCD has been streamlined 
by Shamir |39j . The domain- wall was replaced by Dirichlet boundary conditions. 
The left and right handed components of the Dirac fermion are separated with 
one chirality bound exponentially on one wall, and the other on the opposite 
wall. Using the notations in eqs. (|1I3|I the five-dimensional Dirac operator has 
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the form 

^Dw - \ [iHd'' + d*") - d^d*"] +Dw-mo, (21) 

where Dw is the four-dimensional Dirac operator as defined in eq. ||5J), and 
< Too < 2.^ The Dirichlet boundary conditions set up in 2:5 = and x^ = N 
read 

PrX{xi,X2,X3,X4,,0) = 0, Plx{xi,X2,X3,X4,,N) = 0, (22) 

where Pr/l — (1 i 7^)/2 and x is the five-dimensional fermion field. In the 
limit N —^ 00 the resulting effective four-dimensional theory is chiral symmetric 
on the lattice. I postpone the discussion on this most important point. Let 
me mention that, at the time of its discovery, it was not an issue whether the 
domain-wall idea had anything to do with the Ginsparg- Wilson relation. 

This was a real breakthrough. Although Ddw has an extra dimension (or 
'flavor' space), its structure is quite similar to that of the four-dimensional 
Wilson action which raised the hope for a relatively simple generalization of the 
numerical procedures. The computational cost is increased by a factor 0{N), 
where N should be large, but in problems where chiral symmetry plays an 
essential role there might be no other way to proceed. Concerning the present 
status of applications with domain-wall fermions I refer to the recent paper in 
ref. i4Uj, while for an application with overlap fermions see ref. |41| . 



The fixed-point Dirac operator in QCD solves 
the Ginsparg- Wilson equation 



It was noticed in 1997 that the fixed point Dirac operator in QCD satisfies the 
Ginsparg- Wilson relation |33] . This conclusion follows directly from the saddle- 
point equation eq. (|19|l . Considering the r.h.s. of eq. (|19|l on a very fine lattice 
the fermion propagates over a very fine gauge field configuration U. In this case 
the fermion propagator goes to its continuum limit in this classical equation 
and the only term which is not anticommuting with 7^ is R{V)xB,yB- This 
leads to the equation eq. H13|l (or eq. H12|l 'l with a local R. In addition to being 
classically perfect, the fixed point QCD action has this important quantum 
property. The locality of the fixed-point action follows from general arguments 
of Wilson renormalization group theory^ 

In retrospect it is difficult to understand why this trivial observation was 
not made earlier. The free fixed point Dirac operator was explicitly known for 



^For mo < there exist no massless modes, for mo > 2 doublers appear. 
^Actually, one can optimize the averaging function ui in eq. I19i to make D not only local 
but decaying very fast 1131 . 
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different block transformations 02^ Its physical content is just the statement 
that Dfp is local, the energy-momentum dispersion relation is exactly linear and 
the energy runs in (0, oo) (in the first Brillouin zone) like in a chiral symmetric 
continuum formulation. 

It was also illustrated in 33 that the Ginsparg- Wilson relation implies Ward 
identities from which the standard physical conclusions of chiral symmetry fol- 
low. 



8 Index theorem on the lattice 



It is a simple exercise to show that the solutions of the Ginsparg- Wilson relation 
satisfy the index theorem on the lattice 021 • the zero modes of D are chiral and 
the associated index is a topological invariant which represents the topological 
charge on the lattice. 

For notational simplicity consider the case 2R = 1 in eq. (|12|l and assume 
the hermiticity property in Euclidean space D^ ~ j^Dj^. It is then easy to 
show that the spectrum {A} of D lies on a unit circle around the point z = 1 
in the complex plane with the following properties. The real modes A = and 
A = 2 are chiral (i.e. the corresponding eigenvectors are also eigenvectors of 7^ 
with ±1 eigenvalues). If (jjx is an eigenvector with a complex eigenvalue A then 
(j)x- — "t^4>x and they are orthogonal: ((/)>, 7^ ^a) = 

Define the density 

q{x) = Ur{^'>D{x,x)), (23) 

where the trace is taken in color and Dirac space. Consider now 
Y,q{x) = \it{^'D) = -\Tr[^'{2-D)] 

X 

= -i^(2-A)(0A,7''/'A), (24) 

A 

where Tr denotes trace in all the indices and in the first step we subtracted zero: 
Tr7^ = 0. Using the orthogonality properties discussed above only the A = 
modes contribute giving 

^ q{x) = n_ - n+ , (25) 



^An example was already presented in the Appendix of ref. 1241 refering to the unpublished 
PhD thesis of M. Peskin. In 1982 the community was not yet ripe to appreciate the Ginsparg- 
Wilson paper and pieces of its content were independently rediscovered later. 
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where n_ {n+) is the number of left (right) handed zero modes of the Dirac 
operator. The index of D is, therefore, a sum over x of the pseudoscalar density 
q{x). It can be shown 01] that on smooth configurations q{x) is the continuum 
topological charge density 

9(^) = ^^t^^P<rtT^iF^^Fp^) + O(a') , (26) 

oZtt 

where the trace is taken in color space. In the case of the fixed point Dirac oper- 
ator one can show in addition that ^^ q{x) is the fixed-point topological charge 
Qfp which takes integer values on any gauge configuration, i.e. the topological 
charge from the gauge and from the fermion sector is always consistent on the 
lattice if the FP action is used. 



9 Ginsparg- Wilson fermions from the domain- 
wall construction 



The simplicity of the index theorem demonstrated the power of the Ginsparg- 
Wilson relation. It was obvious that this relation is relevant. In a short time a 
second solution, the overlap Dirac operator, was identified 0S|- Unlike the fixed 
point Dirac operator this solution has a simple explicit structure which allows 
to represent it in a computer to machine precision. In addition, the overlap 
operator binds the Ginsparg- Wilson approach to the domain-wall idea. 

I present here this connection in a form as it occured after several simplifying 
steps |3ni and summarized in ref. |47j. As eq. (|22|l shows the d ~ A fermion 
fields can be naturally identified with 

^.^Plx{x,1)+Prx{x,N), ^^^x{x,1)Pr + x{x,N)Pl, (27) 

where x ^ Z4. The inverse of {ipxi^x) defines an effective four-dimensional Dirac 
operator D^ . In the N —^ 00 limit D^ -^ D can be written as 



aD^l-AiA^'AyK (28) 



where 



A = -a5{Dw - mo) i 1 - -a^iD^ - mo) j (29) 

and a and a^ are the lattice units in the four-dimensional space and in the fifth 
direction, respectively. Taking the limit a^ —>■ one obtains the overlap operator 

D,^ = l-A{A^A)-i, (30) 

where A = — (i^w ^ rno). Both Z?ov and the Dirac operator D in eq. H28|l 
satisfy the Ginsparg- Wilson relation with 2R = 1. This can be shown easily by 
observing that both operators have the form 1 — V, where V is unitary. 



14 



A valid Dirac operator should be local, i.e. it should have an extension of 
0{a) (inverse cutoff). It has been shown 03] that a solution of the Ginsparg- 
Wilson equation can not be ultralocal®. Ultralocality is, however, not necessary. 
Locality (and so, universality) requires that the couplings in the action decay 
with the distance faster than any physical signal. An exponential decay of the 
couplings like exp(— Kr/a), where k is 0(1) defines a local operator. It has been 
shown 05] that the overlap Dirac operator satisfies this bound if the theory on 
the lattice is close to the continuum limit. 



10 Exact chiral symmetry transformation 



As I discussed before, one can raise arguments that the Ginsparg- Wilson relation 
implies the physical consequences of chiral symmetry. One can derive Ward 
identities also and demonstrate that they lead to the same physical predictions as 
those with exact chiral symmetry |50|. This approach is, however, cumbersome. 
In 1998 Liischer has found an exact symmetry transformation which could be 
identified as the chiral symmetry transformation on the lattice [3]. This was 
a theoretical breakthrough. It gave an elegant technical tool to derive Ward 
identities, study anomalies and to handle the problem of chiral gauge theories. 
Chiral symmetry is realized on the lattice in an unusual, but beautiful way. 

Consider a vector gauge theory with Nf flavors and take 2R = 1 in eq. H12|l 
for simplicity. If the Dirac operator D satisfies the Ginsparg- Wilson relation 
then the infinitesimal non-singlet transformation 

i;' = i; + ieTj^ (l- -Dj^j, iJ' = iJ + ie^^ (l - -d] -y'^T (31) 

and the singlet transformation 

iP' = iP + ie-f'^ (l - -d] i; , 11^' =i^ + ielp(l- -Dj-f^ (32) 

leave the fermion action ipDiJj invariant. In eq. 1)31(1 T is an SU{Nf) generator. 
It is trivial to demonstrate this statement. Less trivial is the way this sym- 
metry transformation avoids the Nielsen-Ninomiya no-go theorem: the chiral 
transformation is not a simple 7^ rotation as the no-go theorem assumed. The 
transformation depends on D and so it depends on the gauge field over which 
the fermions propagate. In the formal continuum limit the transformation goes 
over to the standard 7^ rotation^. 

®An ultralocal Dirac operator has a finite number of femiionic offsets. The Wilson Dirac 
operator -Dwi for example, is ultralocal having 9 offsets. 

^In eqs. 1311321 D is multiplied by the lattice unit a as it is obvious from dimensional 
analysis. 
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The action is invariant under the singlet transformation also, but the fermion 
integration measure is not |5H . A non-trivial Jacobi determinant is generated 
depending on the gauge field 

Y[ chpU^''. = (1 + leTvi^'^D)) Yl d^^diP, . (33) 

X X 

For a non-singlet transformation the Jacobi determinant is 1 due to tr(r) = 
0. According to the index theorem on the lattice, a factor of 1 -t- ie2NfQto-p 
is produced by the measure for a singlet transformation, where Qtop is the 
topological charge of the gauge configuration as defined by the index of the Dirac 
operator. Consider the transformation in eq. 1)321) as a change of variable in the 
fermionic path integral. The expectation value (O) of an arbitrary operator 
remains unchanged which leads to the correct anomalous Ward identity on the 
lattice 

{5O)+2Nj{QtopO)=0, (34) 

where SO — [©(i/i ,"0') — ©(V^, '0)]/Je is the variation of the operator O under 
the transformation. 



11 Left and right handed fermions, the mass 
term and the ^ parameter 



In the formal continuum the massless fermion action of a vector theory falls into 
left handed and right handed parts 

'^Dil; = '^lD^^l + l^nDi^R . (35) 

This is also possible on the lattice although in this case the projectors on the tp 
field depend on the gauge configuration |5^I5S1[T^ : 

■0Z, = Phi^ , ipR = Pr^ , V'l = '4^Pr I ■0ij = ^Pl , (36) 
where 

PR/L^lil±l5), PR/L^lil±l5). (37) 

Here 75 depends on the gauge field through the Dirac operator 

75=75(1-^) (38) 

and satisfies (75)^ — 1 due to the Ginsparg- Wilson relation (I consider 2R — 1 
for simplicity). Using these definitions one obtains the decomposition in eq. H35|l 
on the lattice. 

The relations above are striking in many ways. The decomposition ipi^) = 
tpLix) + iPr{x) depends on the gauge field in the neighborhood of x. Even for 
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free fermions, the projectors depend on the Dirac operator. Beyond that, there 
is an asymmetry between the fermions and antifermions: the projectors on ip 
are standard, they are identical to those in the continuum. This asymmetry is 
deeply related to the fermion number anomaly in chiral gauge theories |54| . 

The scalar and pseudo-scalar densities are introduced as usual 

S = IPl^r + ^flV^L = ^ ("l - il? j V , (39) 

P = ^iVii-?fl^L=?75('l-^l?')V'. (40) 

These densities behave correctly under chiral transformations. A mass term in 
the action of a vector theory like QCD can be introduced as 

S = Sg{U) + I dx \^DiP + i^L^tl^R + '^Rm)ipL - ieq{x)\ , (41) 

where a CP breaking 0-term is added also and a sum over flavors is assumed. For 
a real mass m the massive Dirac operator can be written as (1 — m/2){D + m). 
The mass shifts the eigenvalue circle by m and modifies the radius in such a 
way that the point X — 2, where the ultraviolet eigenvalues are concentrated 
remains fixed. Starting from eq. H41(l the known Ward identities can be derived, 
but this time in a controlled environment. The Ward identities are valid at any 
finite value of the cutoff (and so also in the continuum limit) as far as the Dirac 
operator is a local solution of the Ginsparg- Wilson relation. Conserved vector, 
axial vector and chiral currents can be defined jSH EH- It became possible 
to clarify theoretically opaque quantities like the topological susceptibility |58[ 
ICT] and to make progress on the related Witten-Veneziano relation ISU]. The 
Ginsparg- Wilson relation has also been used to construct actions with improved 
chiral and scaling properties, like the chirally improved (CI) action [5 7|. 



12 Chiral gauge theories 



Standard perturbatively defined regularizations have problems with chiral gauge 
theories. Of course, in principle one can give up gauge symmetry and introduce 
all the counter terms needed to get gauge invariant renormalized predictions 
in perturbation theory. This might not be the best way to proceed, not even 
practically. In addition, our ultimate goal is to understand these theories (like 
the standard model) beyond perturbation theory. 

The decomposition eq. (|35|l is the first step towards a gauge invariant for- 
mulation of chiral gauge theories on the lattice. The projectors in eq. H36|l are 
gauge field dependent, so the dividing surface between the left and right handed 
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fermions in the space of fermion degrees of freedom is moving as the gauge con- 
figuration is changing. Even the number of fermions with a certain chirahty, 
say left handed, is changing. Since the number of left handed antifermions is 
constant (their projectors are gauge field independent), the total fermion num- 
ber is changing as we are flying over the gauge configuration space. Indeed, the 
difference between the fermion and antifermion degrees of freedom is 

TtPl - TtPr = ^Tvi-f^D) , (42) 

which, according to the index theorem discussed before, takes different values 
in the different topological sectors. On one hand, eq. 142|l shows that fermion 
number violation occurs naturally in this formulation and this a welcome fea- 
ture. On the other hand, having different degrees of freedom in the different 
topological sectors might imply difficulties to connect these sectors with each 
other in establishing the theory. 

In the last few years important results were obtained in this field. Liischer 
has demonstrated that a U(l) chiral gauge theory with fermions in anomaly 
free representation can be fully defined on the lattice so that gauge symmetry 
is exactly preserved [01]. The same could be established for a general compact 
group in every order of perturbation theory ;fi2. Beyond the f/(l) case it is 
not known whether and under which conditions these theories are free of non- 
perturbative obstructions (anomalies). Such obstructions might exist as has 
been shown in an SU(2) chiral theory with a single left handed fermion in the 
fundamental representation [HJ. This was also demonstrated on the lattice |64|. 
A different approach which is pursued on the lattice since a long time is to give 
up explicit gauge symmetry. I refer here to a recent paper ;65_ and to the 
references therein. 



13 Closing remarks 



The Metropolis algorithm also celebrated its fiftieth anniversary recently |66|. 
Nobody thought fifty years ago that Yang-Mills theory, this beautiful theoretical 
construction, will have so much to do with a stochastic algorithm as it is the case 
since almost three decades. These simulations delivered a lot of non-pcrturbative 
results on the glueball spectrum, static potential and the related fluctuating flux 
tube and thermodynamics. The results are quantitative and some of them are 
quite precise. The numerical results can be connected to analytic predictions 
in some corners of the theory (for a highly non-trivial example I refer to the 
contribution of van Baal in this volume). 

QCD is, of course, more difficult and much more exciting. If the doubts 
concerning the validity of the running staggered fermion simulations |231 are 
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positively clarified we shall certainly see a large number of interesting non- 
perturbative QCD predictions during the next few years. I think, this will 
happen in any case soon. 
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